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Eleven years ago we proposed a model for the electronic configuration of ions of the uranyl type. 
The solution absorption spectra of the (Pu0 2 ) + ^ and (Np0 2 )+ ions are here interpreted on the basis of 
this model. The interpretation is based on exact calculations, and on plausible arguments about the 
strengths and widths of the absorption peaks. If suitable values are chosen for the Coulomb integrals, 
the spin-orbit coupling constant, and the ligand field parameters then the calculated energy levels agree 
well with the observed positions of the absorption peaks. 
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1. Introduction 

Eleven years ago we published a paper [l] 1 on the 
magnetic and spectroscopic properties of the plutonyl 
ion. We now know [2] that the parameters which we 
used in calculating the energy levels were assigned the 
wrong values. Despite this unhappy circumstance 
the correct interpretation was given to the dominant 
peaks in the absorption spectrum and, since the ground 
state of the ion was correctly identified, the calculated 
magnetic properties agreed with the meager available 
experimental data. 

Our purpose in this paper is partly to present revised 
values of the parameters and new calculations of the 
energy levels. Our primary purpose is to discuss in 
detail the absorption spectra of the isoelectronic 
(Pu0 2 ) ++ and (Np0 2 )+ ions, both of which have two 
electrons outside what we might call the uranyl core. 
We have been curious to see the extent to which plaus- 
ible, semiquantitative arguments about strengths and 
widths of absorption lines could be used in identifying 
the rather complex solution absorption spectra of 
these ions. 

In analyzing the (Np0 2 ) + and (Pu0 2 ) ++ spectra the 
first important question which one has to decide is 
whether the ligand field strengths for these ions are 
the same as those for the neptunyl ion. Information 
bearing indirectly on this point comes from the infrared 
measurements of Jones and Penneman [3] on uranyl 
and transuranium (V) and (VI) ions in aqueous solu- 

1 Figures in brackets indicate the literature references and notes at the end of this paper. 



tions. Table 1 gives the asymmetric stretching fre- 
quencies of these ions as determined by Jones and 
Penneman. The similarity of the stretching fre- 
quencies for the first four ions in the table indicates 
that the ligand fields are similar for all four. On the 
other hand, since for (Np0 2 ) + and (Am0 2 )+ the vibra- 
tion energies are lower, the bond length is probably 
greater and the ligand field weaker. 

TABLE 1. Asymmetric stretching frequencies of uranyl-type ions 
in aqueous perchloric acid solution 



Ion 


Frequency/c 




cm " ' 


(UOt) ++ 


965 


(N P 2 ) + + 


969 


(Pu0 2 ) ++ 


962 


(Am(),)+ + 


939 


(Np0 2 ) + 


824 


(Am0 2 ) + 


832 



Indirect information on the variation of the ligand 
field strengths in the uranyl type ions can also be 
obtained by studying the variation of the crystal field 
parameters in any series of structurally similar com- 
pounds. For example, for several rare earth trichlo- 
rides [4] the variation from one rare earth ion to the 
next is probably 10 or 20 percent. 

Data on the absorption spectrum of the (U0 2 ) + ion 
would yield important information on the strength 
of the ligand field in (Np0 2 ) + if it is really the case that 
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similar ions have similar fields. The only experiments 
on the (U02) + absorption have been done by Kraus, 
Nelson, and Johnson [5]. They found no absorption 
peak with a molar extinction coefficient as high as 1.5 
in the region between 10,000 and 27,800 cm -1 . The 
peaks corresponding to the transitions to |1 — > and 
1 1 + > should fall in this region. In neptunyl they 
have extinction coefficients of approximately 6. The 
more intense absorption peaks, which the (UO2K 
spectrum should have, occur below 10, 000 cm -1 . 

To begin the calculations we assumed that the ligand 
field strengths were the same for (Pu02) ++ as for 
(NpC>2) ++ and were about 15 percent smaller for 
(Np0 2 ) + . This assumption proved to be a reasonable 
starting point. 

One expects the value of £ to be larger for the 
plutonyl ion than for the neptunyl ion because the 
spin-orbit coupling constant generally increases with 
increasing atomic number. £ may not have quite the 
same value for (Np0 2 ) + as for (Np0 2 ) ++ because the 
overlap of the wave functions on the surrounding ions 
may be slightly different. 

The important additional complication in the two- 
electron spectra is the Coulomb repulsion of the 
electrons. The calculation of the matrix elements 
is considerably facilitated by use of table 1 of 
Condon and Shortley's [6] book. All the matrix ele- 
ments of the Coulomb interaction depend on numerical 
factors, which can be obtained from the table, and on 
the four radial integrals F (0) , F (2) , F {4 \ and F (6) or, 
more conveniently, the parameters Fo = F (0 \ F> 
= F< 2 7225, F 4 = F (4) /1089 and F« = 25 F< 6 >/(429) 2 . 
In practice, we have followed Racah by using the linear 
combinations: 

/. = 5F 2 + 6F 4 -91 F e 
A = 924<F« 
B=U F 4 -MF H . 



We used a Fermi-Thomas model to calculate a 5/ 
radial wave function appropriate to Pu 5+ from which 
Basu calculated F (2) = 79,000 cm -1 . We then used an 
approximate wave function to calculate the relative 
values of the F (2w) . We found, for example, that 
L = 1577 cm -1 . A further useful guide comes from 
the wave functions calculated by Ridley [7] for U 6+ by a 
Hartree method. For the 5/ electrons she found 
F = 137,340 cm- 1 , F 2 = 291.9 cm" 1 , F 4 = 39.68 cm 1 
and F 6 = 4.281 cm -1 ; therefore L= 1338 cm -1 . These 
calculations give one only a rough idea of the values of 
L, A, and B. The actual values of these parameters in 
(Pu0 2 )+ + and (Np02) + must be determined by analyzing 
the absorption spectra. 



2. Procedure 

In our analysis of the (Np02) ++ data in III we had ten 
parameters at our disposal when covalent effects and 



ligand fields of threefold and sixfold symmetry were 
considered. For (Np0 2 ) + and (Pu0 2 ) ++ six additional 
parameters are required to characterize the Coulomb 
interaction (covalent effects included). The total of 
16 parameters is inconveniently large. We have 
therefore simplified the problem by leaving out of con- 
sideration covalent effects and ligand fields with three- 
fold symmetry. We hoped that by varying the eight 
remaining parameters we would be able to obtain 
agreement with at least the gross features of the ob- 
served spectra, and this hope has been fulfilled. It is 
not difficult to include the neglected effects in the cal- 
culations, and this should probably be done when 
low-temperature absorption spectra of crystalline 
uranyl-type compounds become available. 

In II we excluded from consideration the cr-states — 
i.e., states in which mi for one or both electrons is 
zero. However, we have included these states here. 
When the Coulomb and the spin-orbit interactions 
are taken into account, and the ligand fields of axial 
and sixfold symmetry, one obtains secular equations of 
order 12, 5, 14, 16, 7, and 7. These numbers corre- 
spond, respectively, to crystal field quantum numbers 
0, 0, 1, 2, 3, and 3. The secular equations can readily 
be solved, for particular values of the parameters, by 
using an electronic computer. 

Our procedure in analyzing each spectrum was to 
select values for all the eight parameters and then to 
calculate eigenvalues and eigenvectors. The cal- 
culated eigenvalues were studied in relation to the 
known absorption spectrum and an attempt was made 
to find improved values of the parameters. These 
improved values were then used for another calculation 
of the energy levels. The procedure was repeated 
until we had what we felt was the most satisfactory 
possible interpretation of the spectrum. 



3. The Plutonyl Spectrum 

For both the (Np0 2 ) + and the (Pu0 2 ) ++ ions the 
principal component of the ground state is the 3 H 4 
1 3 — , 2 — > state [8, 9]. In the absence of the "crystal- 
line" field Hund's rules place this state lowest. In 
all the calculations we have done, with widely varying 
values of the parameters, 3 H 4 has been found to be 
lowest. Finally, and most important, the paramag- 
netic resonance data on crystalline plutonyl rubidium 
nitrate [10] and plutonyl sodium acetate [11] indicate 
unambiguously that the ground state is 3 H 4 . There- 
fore in discussing the absorption spectra we shall 
assume, without further comment, that all transitions 
take place from a ground state which is predominantly 
3 H 4 . 

The absorption spectra of the (Pu02) ++ and (Np0 2 ) + 
ions are shown in figures 1 and 2. It will be seen that 
the two spectra have certain similarities particularly 
in their low frequency regions. The spectrum of 
(Pu02) ++ is shifted towards higher frequencies relative 
to that of (Np0 2 ) + . 

The most prominent peak in the plutonyl spectrum 
is the one at 12,037 cm -1 for which the molar extinction 
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coefficient is 300, the width at half maximum is 34 cm -1 
and the strength [12] is about 31. We regard the 
identification of this peak as predominantly 3 n 2 
|2 — , —3— > as certain. The transition which takes 
place is |3-, 2->-»|-3-, 2->. We have al- 
ready argued in III that the 3— > — 3 transitions should 
give rise to narrow absorption peaks. The fact that 
they are strong as well was of crucial importance in 
the interpretation of the neptunyl spectrum. By 
suitably choosing the values of the Coulomb integrals 
and the spin-orbit coupling constant other levels, 
rather than the |2 — , — 3 — > level, can be made to 
fall at 12,037 cm -1 . However, if such a choice is made 
one is left without a plausible identification for the 
next most prominent peaks at 10,185 and 10,500 



1 

25 — 
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15 — 

1 
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We identify these peaks, in part, as due to the transi- 
tions 2^ — 2 to a split 3 n level. The apparent 
degeneracy of the states |3 — , — 2 — > and |2-f, 
— 3+> is removed by the spin-orbit coupling in a 
second order perturbation calculation [13]. The 
symmetric combinations 



13- 



- 2 - > + | 2 + , - 3 + > 



Figure 2. Absorption spectrum o/Np (V) in 1.0 M DC10 4 — D 2 0. 

The curve has been replotted from data given by W. C. Waggener in J. Phys. Chem. 
62,382 (1958). 



etc., are coupled by £ L • S to the 2 states. The 
antisymmetric combinations, however, are not coupled 
to any other states by £ L • S. They are coupled to 
one another by the Coulomb interaction, and to 1 Iq 
by the ligand field component which has sixfold sym- 
metry. We shall refer to these antisymmetric 
combinations as 3 Ho states, and to the symmetric 
combinations as 3 IIo states although this notation is 
the reverse of that used in our previous paper on the 
plutonyl ion (II). One finds, as a result of calculation, 
that the 3 IIo |3 — , — 2— > state lies below the corre- 
sponding : Tlo state. Since 3 II ( , ) |3 — , — 2 — > interacts 
weakly with the other 3 IIi states, and since it can be 
reached from the 3 H. 4 ground state by a 2— » — 2 transi- 
tion which leaves the charge distribution nearly un- 
affected, we conclude that the ;i II peak should be 
narrow. If we adjust the parameters of our theory so 
that 3 I1 2 is located near 12,000 cm -1 then we find that 
: *n ( ') lies near 10,000 cm -1 . Its position is rather 
insensitive to all parameters except £ and L. We 
therefore identify the peak at 10,185 cm -1 as due to 
transitions to 3 Ili, and we regard this identification 
as virtually certain. 

The experimental peak at 10,500 cm -1 probably 
consists of two unresolved peaks since the absorption 
curve has a well developed shoulder [14] at about 
10,700 cm -1 . One of these unresolved peaks is due 
to transitions to 3 I1() |3 — , — 2 — >. The peak should 
be wider and also less intense than the 3 I1 peak 
because of the extensive coupling of 3 Ilo to other states. 




FIGURE 1. Absorption spectrum of Pu (VI) in 0.5 M HC1. 

The curve has been replotted from data given by R. E. Connick, M. Kasha, W. H. McVey 
and G. E. Sheline in The Transuranium Elements, National Nuclear Energy Series, IV, 
14b, McGraw-Hill Book Co., New York, 1949. (Paper 4.20, figs. 12-15.) 
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The calculations place 3 n o in the vicinity of 10,700 
cm -1 , and since the shoulder has less strength than 
the peak we have identified with 3 IIo, we ascribe the 
unresolved 10,700 cm -1 peak to 3 n o . 

The other component of the unresolved doublet 
is probably a Si state. The transitions that give rise 
to this absorption peak are |3 — , 2 — >— >|— 2— , 
2-> and |3-, 2->-» |3-, -3->; that is, they 
are 3 — > — 2 and 2-^ — 3 transitions. For the neptunyl 
ion the 3-^ — 2 transitions are almost as narrow 
as the 3-^ — 3 transitions, so one can expect the same 
to be true for the plutonyl ion. They are, however, 
less probable [15]. It therefore seems likely that the 
experimental peak at 10,500 cm -1 should be identified 
with transitions to Si. The principal difficulty with 
this assignment is that the peak has at least twice the 
strength we would expect from our study of the 
neptunyl spectrum. We regard the identification of 
the unresolved structure between 10,300 and 10,800 
cm -1 as partly Si and partly 3 n as pretty certain, 
but it is difficult to say precisely where the two peaks 
should lie. 

The two small peaks at 12,660 and 12,890 cm -1 may 
be due either to the 3 H« and 3 H^ states or to the vibra- 
tional fine structure of the 12,037 cm -1 peak. The 
charge distribution in the 3 Hs states is such that the 
corresponding absorption peaks should be narrow. 
They should also be weak since the absorption is due 
principally to the contamination of the 3 H« states by 
3 IIo states. We consider it unlikely that the two peaks 
are due to transitions to the 3 H 6 and 3 H/; states because 
the separation of these states is less than 100 cm -1 
in almost all of our calculations. 

If the peaks are to be ascribed to vibrations the vibra- 
tional energy would have to be 620 or 850 cm -1 . 
The energies of the symmetric and antisymmetric 
stretching vibrations in the ground state of solid sodium 
plutonyl acetate are 818 and 930 cm -1 respectively 
[16]. Therefore we consider it reasonable to assign 
the peak at 12,890 cm -1 to the symmetric stretching 
vibration. Since the charge distribution in the excited 
3 n 2 state is nearly the same as in the ground state the 
vibrational energies should be nearly the same in the 
excited state as in the ground state. Despite the fact 
that the 3-^ — 3 transitions do not excite many vibra- 
tions the peak at 12,890 cm -1 shows up because of 
the very great strength of the 12,037 cm -1 peak. 

The 3 Hk doublet is probably not resolved, both com- 
ponents being concealed in the single hump at 12,660 
cm -1 . One can perhaps understand in this way why 
the 12,660 cm -1 peak is broader than one would expect 
for 3 H«. 

There now follows a discussion of the identification 
of the remaining peaks of the experimental absorption 
spectrum of the plutonyl ion. These identifications 
rest on arguments about intensities and breadths and 
on detailed calculations of the positions of the energy 
levels. In contrast to the peaks just discussed we 
do not regard the identification of any of the remaining 
peaks as certain. The difficulty we have encountered 
can be stated in the following way. Depending on the 



choice of parameters either there are too few calculated 
levels in the region from 15,000 to 25,000 cm -1 , or 
else there are calculated levels with no corresponding 
observed levels. 

This circumstance makes it necessary for us to 
concentrate our attention on what we think are the 
principal transitions, and on the dominant features of 
the experimental spectrum. Consider the possible 
"allowed" transitions from the 3 H 4 |3 — , 2 — > ground 
state. They are 3^ — 3, —2, — 1, 0, or 1 and 2 
— * — 3, — 2, — 1, 0, or 1. We have already located the 
3^> — 3 and the 2-^ — 2 transitions. The 3-^0 and 
2-^0 transitions are almost certainly well above 25,000 
cm -1 since an electron in an mj=0 state is strongly 
repelled by the axial charge distribution. Due to the 
extensive mixing of the II i and Si states it is unlikely 
that the 3— > — 2 and 2-^ — 3 transitions will give rise 
to any peaks that are stronger than the peak at 10,500 
cm -1 which we have associated with a Si state. The 
remaining possible transitions are 3— >— 1 or 1 and 
2^—1 or 1. The calculations place the excited 
states in the following order of ascending energy: 
3 T 3 and 3 n (2-»l) [17], 3 <D 2 (3->l), 3 A, (2-*-l), 

3 n ( ; and 3 n () (3^-i). 

For the neptunyl ion the 3^2 transition is relatively 
narrow, probably because the change in the charge 
distribution is mainly in the equatorial region and does 
not excite vibrations of the ion, or even of its surround- 
ings, to any large extent. However the 3— »1 and 
3— >— 1 transitions are broad because vibrations are 
excited. It follows that the redistribution of charge 
which accompanies a 2^1 transition also ought to 
excite vibrations. Therefore the absorption peaks 
corresponding to transitions to the 3 r ;i , ; T' } , 3 4> 2 , 
3 Ai, 3 rT ( j, and 3 n states of the plutonyl ion should all 
be broad. By analogy with the neptunyl ion the widths 
at half maximum should be about 1000 cm -1 [18]. 

The first prominent feature of the absorption 
spectrum above 14,000 cm -1 is the sharp peak at 
16,075 cm -1 . Because of its sharpness this peak 
cannot be associated with any of the states listed in 
the previous paragraph. It is a moot question, but 
one which is crucial for our analysis, whether or not 
the 16,075 cm -1 peak rests on a broad base. If so, 
the strength of the underlying peak may be as large 
as 5. It would then be reasonable to associate this 
broad peak with transitions to 3 r :} and 3 r.' } . 

A serious difficulty then arises: It becomes neces- 
sary to associate the strong, broad peaks at 19,100 and 
21,200 cm -1 with transitions to 3 <f>2 and 3 Ai, respec- 
tively. The observed separation of the peaks is 2100 
cm -1 but the calculated separation is about 4000 cm -1 
and cannot be significantly reduced by any reasonable 
adjustment of the parameters. 

Because of this difficulty [19] we consider it more 
likely that the nearly degenerate 3 r :} and 3 T^ states 
should be identified with the broad peak at 19,000 
cm -1 . The principal difficulty now is that there are 
only three calculated levels in the region between 
14,000 and 19,000 cm -1 but the absorption curve indi- 
cates the presence of transitions to at least four states. 
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TABLE 2. Observed and calculated absorption spectrum of th 


e(PuO,) ++ ion 








Position in cm ' 








Composition of wave function 
















Width at half 




hy \ni(\ value 






Level 


Purity" 






^max. 


max. in cm ' 


Strength b 


Calculated Btrength c 








Calc. 


Exptl. 








3 


2 


1 









:, H, 


0.9249 














0.9526 


1.0180 


0.0295 


0.0000 






So 


.7331 


2445 






















Hi 


.5929 


4258 






















'H, 


.9825 


7133 






















:, n;<3, -2) 


.9580 


10157 


10185 


10.2 


150 


- 5.7 


.9617 


0.9993 


.0385 


.0005 


0.8860 (2 -* - 2) + 0.0350 (3-* - 1) 




X, 


.5966 


10489 


10500 


22.2 


195 


16. 


.6953 


1.2645 


.0395 


.0008 


.4091 (3 -*-2) + .1499 (2-* -3) +.0216 (2-*- 


-1) 


:, n<> (3, - 2) 


.628 1 


10610 


10700 


-4.4 


~ 160 


-2.5 


.6518 


1.3295 


.01 1.; 


.0043 


.5916(2-^-2) 




:! 1I, (3, -2) 


.9204 


11892 


1 2037 


300. 


34 


31. 


.9579 


0.9623 


.0788 


.0009 


.85 17 (3 — -3) 




:, H« 


.9278 


12862 


12660 


2.7 


400 


. 


1.0680 


.9318 


.0000 


.0000 


<0.01 




:, h; 


.9580 


1 2866 


3. 


1.0686 


.9312 


.0002 


.0000 


<0.01 




T, (2, 2) 


.8068 


15469 


15420 


0.6 


~ 200? 


- 0.3? 


0.1130 


1.7539 


.1330 


.0000 


.0931 (gnd) 




So 


.9720 


15980 


16075 


~ 5. 


- 135 


- 1.6 


1.1555 


0.7877 


.0564 


.0004 


.0229 (3 -*-l) 




1, 


.6328 


17877 


17800 


2.7 


-600 


- 3. 


1.2731 


.6989 


.0267 


.0012 


.1540(3-* -2) + .4301 (2-* -3) 




3 r 8 


.8828 
.8838 


19080 
19086 


19100 


11. 


~ 1100 


-24. 


0.8978 
.8981 


.1022 
.1018 


.9840 
.9842 


.0159 
.0157 


.8170(2-* 1) 
.8179(2-* 1) 




\ (l 


.8102 


19780 


19810 


~7. 


-80 


- 1. 


.8681 


1.0141 


.1136 


.0042 


.0901 (3^-1) + .0999 (2-* -2) 




3 * 2 (2, 1) 


.8584 


2 1 253 


21200 


14. 


- 1200 


-30. 


.0959 


0.8874 


.9957 


.0209 


.7938(3-* 1) 




'H, 


.8935 


21885 


21840 


~ 7. ? 


-300 


1. 


1 .0028 


.8986 


.0985 


.0000 


.0120(1-* 2) + .0407 (2-* 3) 




Hi 


.6660 


23063 










1.1679 


.7183 


.1125 


.0011 


.0227 (3 -* - 2) + .2926 (2-* - 3) 




A 2 


.7349 


23153 










0.7780 


.2504 


.9580 


.0136 


.1037(3-* J) + .0380 (3-* -3) 




:, r, 


.838 1 


24527 










.8715 


.1750 


.9453 


.0052 


.0229 (gnd) 




■u 


.9301 


24780 










1.9294 


.0704 


.0001 


.0000 


.0125(1-* 3) 




■i; 


.9310 


24781 










1.9310 


.0689 


.0000 


.0000 


.0125 (-*> 3) 




3 A, (3,-1) 


.8506 


24917 










0.9302 


.1416 


.9233 


.0048 


.7870(2-*-!) 





;| Hy "purity" we mean the extent to which the level can he characterized hy the designation in colun 
b The formula used in computing the strength is given in footnote 12. 
' Sec the appendix. 



When the 3 r ;} states are placed at 19,000 cm -1 the 
broad peak at 21,200 cm -1 can readily be identified with 
transitions to the 3 <J> 2 |2 — , 1 — > state. 

With the dominant peaks in the spectrum at least 
tentatively identified we can refine out estimates of the 
various parameters, repeat the calculations, and then 
try to identify the remaining absorption peaks. 

It appears that the narrow peaks at 16,075 and 19,810 
cm -1 should be associated with transitions to states 
that are predominantly So. This assignment is reason- 
able in view of the positions, widths and strengths of 
the two peaks. 

We associate the peaks at 15,420, 17,800, and 21,840 
cm -1 with transitions to T4, Oi, and l H5 states, respec- 
tively. These assignments are all dubious and are 
made only because there are no alternatives. l T 4 is in 
the right position but it appears fairly broad in figure 1 
whereas the calculations indicate it should be narrow. 
Also its calculated strength would be negligible were it 
n6t for the fact that its wave function contains some of 
the same zero-order states as the ground state wave 
function. 

The Si state whose calculated position is 17,877 
cm -1 should be narrow and perhaps half as strong as 
Si at 10,500 cm -1 . It is possible that the 111 state 
is an unresolved part of the large absorption peak at 
19,100 cm- 1 . 

The difficulty with the *H5 assignment is that the 
calculated strength seems to be too small. 

We have now discussed all calculated levels from 
10,000 to 23,000 cm -1 . There is however one more 
observed peak at 16,950 cm -1 . We are not entirely 
convinced that this peak exists since it seems to be 
located by a single point on the most detailed figure 
in the paper of Connick et al. [20], and disappears 
when the HC1 concentration is increased. 



It is possible that both the 16,950 (if real) and the 
17,800 cm -1 peaks are associated with the major peak 
at 19,100 cm -1 in the same way as the ripples on the 
low energy side of the neptunyl absorption peaks are 
associated with the strong 19,100 and 21,100 cm -1 
peaks. 

The observations and the calculations on the plutonyl 
ion are summarized in table 2. See the appendix for 
an explanation of the right hand side of the table. 



4. The (Np0 2 ) + Spectrum 

The (Np0 2 ) + ion, like the (NpO,) ++ ion, is believed 
to be linear. The optical absorption spectrum of the 
(Np0 2 ) + ion in aqueous solution has striking similari- 
ties to that of the (Pu0 2 ) ++ ion so we shall assume that 
(Np0 2 ) + is an ion of the uranyl type and attempt to 
interpret its properties on this basis. It has already 
been mentioned that the asymmetric stretching 
frequency for the (Np0 2 ) + ion is lower than for the 
neptunyl or plutonyl ions. Also, the (Np0 2 ) + spec- 
trum, as compared with the plutonyl spectrum, is 
shifted toward the red. We attribute these shifts to a 
weakening of the molecular field, to a lower value of the 
spin-orbit coupling constant, and to lower values of 
the Coulomb integrals for the (Np0 2 ) + ion. 

The solution absorption spectrum of the (Np02) + ion 
is shown in figure 2. We regard the identification 
of the following absorption peaks as nearly certain: 

8953 cm- 1 3 n o 

9146 Si 

9780 3 n 

10208 3 n 2 
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TABLE 3. Observed and calculated absorption spectrum of th 


e (Np0 2 ) + ion 






Position 


in cm -1 








Composition of 


wave function 














Width at half 






by \mi 


value 






Level 


Purity" 






e max 


max. in cm -1 


Strength 1 ' 










Calculated strength 1 




Calc. 


Exptl. 






3 


2 


1 







:, H 4 


0.9216 














0.9626 


0.9991 


0.0382 


0.0000 




2o 


.6920 


2475 




















n. 


.6181 


3893 




















:i H 5 


.9115 


6148 


6173 


3.5 


47 


1.2 


.9692 


1.0093 


.0215 


.0000 


.0333 (3 ->- 2) + 0.0186 (1-* 2) +0.8400 (+-*-) 


3 n^(3,-2) 


.8380 


8725 


8953 


5.8 


55 


1.3 


.9483 


.9943 


.0570 


.0003 


.7723 (2->-2) +.0523 (3-»-l) 


:! X, 


.6492 


9450 


9146 


22.4 


137 


13. 


.8301 


1.1163 


.0530 


.0006 


.3704 (3-* -2) + .2332 (2-> -3) + .0350 (2-> -1) 


3 n (3,-2) 


.4804 


9862 


9780 


7.7 


83 


2.5 


.7417 


1.2427 


.0129 


.0027 


.4486 (2-* -2) 


:i n 2 (3, -2) 


.8503 


10294 


10208 


>82 


< 139 


45? 


.9250 


.9311 


.1427 


.0011 


.7851 (3->-3) 


vib 






10970 


2.1 


93 


0.7 












:i He 


.8225 


11171 


11160 


0.5 


~ 150 


-0.3 


1.0627 


.9274 


.0098 


.0000 


.0876(2-^-2) 


:{ H 6 


.7076 


11199 


1.0213 


.9753 


.0028 


.0006 


.1352 (2-* -2) 


:i r.; 


.9122 
.9162 


12949 
12969 


13020 


2.1 


-900 


-9. 


.9301 
.9315 


.0700 
.0685 


.9901 
.9908 


.0099 
.0091 


.8410(2^ 1) 
.8447 (2-*- 1) 


lo 


.9006 


13894 


13824 


<0.5 


-200 


<0.3 


1.1547 


.6680 


.1754 


.0020 


.0739 (3->-l)+.0103(l->-3) 


T 4 (2, 2) 


.6550 


14345 


14577 


2.4 


80 


0.5 


.1897 


1.5267 


.2835 


.0001 


.0814 (gnd) 


:i S, 


.5931 


15911 


16220 


~-9 


160 


>5.5 


1.0743 


.8651 


.0592 


.0014 


.2132 (3-> -2) + .3155 (2-> -3) 


s *« (2. 1) 


.6598 


16045 


16100 


<7.6 


-700 


- 11. 


.2906 


.7095 


.9838 


.0159 


.6081 (3-* 1) + . 0232 (3-> -3) 


A, (3,-1) 


.4482 


16907 


16906 


-1.6 


-150 


-0.5 


.6535 


.4065 


.9355 


.0047 


.2733 (3-* 1)+ .0451 (3->-3) 


So 


.7912 


17169 










.7888 


1.0511 


.1570 


.0030 


.1288 (3-»- 1) + .0524 (2-> -2) 


:i Ai 


.5633 


18085 


18116 


2.3 


- 1000 


-5. 


.9953 


.2706 


.7339 


.0003 


.0318 (2->-3) + .5202 (2-* -1) 


3 r 4 


.6659 


18149 










.8568 


.2518 


.8866 


.0047 


.0229 (gnd) + .0125 (3-* -3) + .0135 (3-> 1) 


'H 5 


.3654 


18589 










.9410 


.5942 


.4620 


.0029 


.0585 (2-* -3) + .2540 (2-* - 1) + .0128 (2-> 3) 


n, 


.5895 


19452 


19360? 


>0.3 


-500 


-0.3 


.7486 


.7760 


.4707 


.0047 


.0902 (2-> -3) + .0271 (2-» - 1) 


:i <P :t (2. 1) 


.8810 


20337 










.0993 


.9008 


.9937 


.0064 


.0394(2-^ l)+.0182(3->2) 


3 <J>3 (2, 1) 


.9240 


20425 










.0554 


.9446 


.9961 


.0039 


.0509 (2-* 1)+ .0191 (3-* 2) 


'Ie 


.9216 


20532 










1.9205 


.0785 


.0012 


.0000 


.0140(1-* 3) 


T 6 


.8782 


20578 


21008 


6.4 


320 


3.7 


1 .8369 


.1580 


.0049 


.0003 


.0133 (l->3) 


n, 


.7665 


21359 


.5242 


.8574 


.6020 


.0164 


.1465 (2-* -3) + .0498 (3-> -2) 


:i A.;(2, i) 


.9811 


21841 










.9976 


.0024 


.9907 


.0093 


< .01 


:! A :i (2. 1) 


.9289 


21912 


21700? 


3.5? 


-800 


-5. 


.9525 


.0475 


.9925 


.0075 


.0155 (2-> 1) 


3 no(2,-i) 


.9220 


22380 


22600 


5.6 


-700 


-7. 


.0664 


.9872 


.9342 


.0122 


.8497 (3->-l) + .0595 (2 -*-2) 


:! n (2, -l) 


.6347 


23397 


23400 


8.7 


-800 


-11. 


.2989 


.9586 


.7327 


.0098 


.5858 (3-»- 1) + .0582 (2-> -2) 


T 5 


.7056 


24292 










.9954 


.2835 


.7212 


.0000 


.0102 (2->-l) 


:, n, (2. -l) 


.7623 


24399 










.2109 


.7964 


.9779 


.0149 


.0191 (3-»-3) 



a Bv "purity" we mean the extent to which the level can be characterized by the designation in column 1. 
''The formula used in computing the strength is given in footnote 12. 
c See the appendix. 



The reasons for these assignments are the same as 
those for the corresponding assignments for the plu- 
tonyl ion. There are two small difficulties here which 
keep us from being completely certain about the 
assignments. (1) The peak at 9146 cm -1 has a curious 
top [2] which makes it look as though it might be two 
unresolved peaks. (2) The 3 n peak is stronger than 
the 3 n o peak. Since 3 I1 is more nearly a pure state 
the reverse should be true, as is the case for the 
plutonyl ion. However, it is difficult (impossible ?) to 
find an interpretation of the peak which is alternative 
to the one given above. 

The peak at 6173 cm -1 is due primarily to the mag- 
netic dipole transition from 3 H 4 to 3 H 5 . The strength 
of this peak, calculated according to the formula given 
in III is 1.2; for a pure magnetic dipole transition the 
strength should be 1.0. It is rather unusual to find 
a line in a solution absorption spectrum that can be 
identified as due to a magnetic dipole transition. The 
calculated position of the 3 H 5 peak is very near the 
observed position so we regard the identification of 
this peak as certain. 

There are two small peaks at 10,970 and 11,160 
cm -1 in the (Np0 2 ) + spectrum. If one of them is 
assigned to the 3 H 6 states and the other to excitation 
of a vibrational mode the energy associated with the 
vibration would have to be 760 or 950 cm -1 . We do 
not have any direct information on the symmetric 
stretching frequency for the (Np0 2 ) + ion in aqueous 
solution. However we note that according to the data 



of Jones [16], and of Jones and Penneman [3], the ratio 
of the asymmetrical stretching frequencies for (Np0 2 ) + 
and (Np0 2 )++ in solution is 824/969 while the ratio 
of the asymmetrical stretching frequencies in solution 
and in a crystalline solid is 969/934. Therefore we 
estimate that the energy of the stretching vibration 
should be (824/969) X (969/934) X (energy of symmetric 
stretching vibration in a crystalline solid). The latter 
number [16] is 844 cm -1 so our estimate for the energy 
of the symmetric stretching vibration of the (Np0 2 ) + 
ion in aqueous solution is 745 cm -1 . We therefore 
associate the peak at 10, 970 cm -1 with a combined 
electronic and vibrational transition and the peak at 
11, 160 cm -1 with the transitions to the 3 H 6 and 3 H' 6 
states. 

As in the case of the plutonyl ion we do not regard 
our identification of any of the remaining peaks in 
the (Np0 2 ) + spectrum as certain. Generally speak- 
ing, our procedure in this case was the same as before: 
we interpreted the relatively strong, broad peaks as 
due to the "allowed" transitions to the 3 r 3 , 3 <J> 2 , 3 Ai, 
and 3 n (2 — 1) states. We found that parameters 
could be chosen so that the calculated positions of 
these states agreed with the observed positions of 
strong, broad peaks insofar as they could be 
determined. 

We believe one of these peaks definitely exists and 
is centered at 13, 020 cm -1 . We ascribe it to transi- 
tions to the 3 r 3 and 3 r 3 states. We place the next 
broad peak somewhat arbitrarily at 16, 100 cm -1 . It 
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lies under the fairly sharp, strong peak at 16, 220 
cm -1 . We associate the broad peak with the 3 ^>2 
state. The plots we have made with expanded ab- 
scissa reveal the structure much more clearly than 
figure 2 does. 

The third broad peak is at about 18, 120 cm -1 . We 
identify it as being due to transitions to the 3 Ai state. 

The peak at 21,010 cm -1 is too narrow to be asso- 
ciated with a major redistribution of charge. It ap- 
pears to rest on the low energy side of a broad peak 
centered at about 21, 700 cm -1 . There are two more 
broad peaks at about 22, 600 and 23, 400 cm- 1 . We 
identify two of these three peaks with transitions to 
the 3 n^ (2, -1) and 3 n (2, -1) states. The most plau- 
sible identification of the remaining peak is achieved 
if the 3 I1() and 3 n states correspond to the peaks at 
22, 600 and 23, 400 cm -1 , respectively. The remaining 
peak at 21,700 cm -1 must then be ascribed to transi- 
tions to 3 A 3 and 3 A,3 states. The principal objection 
to these assignments is that the experimental strength 
of the 21, 700 cm -1 peak is greater than the calculated 
strength for transitions to 3 A 3 states. 

Another possibility is that some of the structure 
observed above 20, 000 cm -1 is due to excitation of 
the vibrational modes of the (Np02) + ion. 

With all the broad peaks now identified plausibly 
but not certainly we consider the remaining features 
of the (Np02) + absorption spectrum. We ascribe 
the moderately strong 16, 220 cm -1 peak to transitions 
to a Si state and the 21,010 cm -1 peak to transitions 
to a 111 state. These assignments are reasonable in 
terms of calculated position, strength and charge 
distribution, and in relation to the width and strength 
of the 3 2i peak at 9146 cm -1 . 

We identify the peak at 14,577 cm" 1 with the T 4 
(2, 2) state. The corresponding transition has a cal- 
culated strength different from zero only because the 
T4 state contains the same zero order wave functions 
as the ground state. The 14,577 cm -1 peak is nar- 
rower than we expect for T4 on the basis of the cal- 
culated charge distribution. The identification is 
therefore, rather unsatisfactory. However, the only 
other state in the vicinity is a So. It is difficult to 
adjust the parameters so that So lies in the correct 
position. 

We ascribe the two small peaks at 13,824 and 16,906 
cm -1 to transitions to So and A 2 states, respectively. 
Again, the assignments are reasonable in terms of 
calculated position, strength and charge distribution. 

The very small peak at 19,360 cm -1 may be due to 
transitions to the 'Hs state or a ITi state. 

5. Values of the Parameters 

The values of the parameters used in the final cal- 
culations on (NpOo) + and (Pu02) + + are given in table 4. 
It is gratifying that X, Y, Z, and V 6 have values that 
are consistent with expectations based on our analysis 
of the (Np0 2 )+ + spectrum in III. Also, the Coulomb 
integrals are larger for the (PuOi) ++ than for the 
(NpOi) + ion. This behavior is what one would expect 
since the higher nuclear charge for (PuOi) ++ pulls 



the electrons in toward the center more and therefore 
increases their Coulomb repulsion. 

The parameters £ and L are probably determined to 
within 5 percent. The levels below 13,000 cm -1 that 
we regard as almost certainly identified depend prin- 
cipally on £ and L and enable us to fix their values 
fairly closely. A is probably determined to within 
10 percent by minimizing the mean square deviation 
for as many levels as we can reasonably identify. 

Table 4. Values of the Parameters 





(Np0 2 ) + 


(Pu0 2 ) + + 




cm~ x 


cm' 1 


X 


253 


- 1183 


Y 


12034 


16793 


Z 


48555 


49650 


V, 


807 


168 


£ 


2090 


2467 


L 


919 


934 


A 


4322 


5926 


B 


184 


36 


F, 


219.5 


254.2 


F 4 


41.21 


41.05 


F 6 


4.677 


6.413 



The value of B is determined largely by the position 
of the i r 4 (2,2) level. As indicated above we are not 
at all sure we have located this level. We therefore 
estimate that B is known only within a factor of 2. 

We started the calculations with V§ ~ 1000 cm -1 . 
For both (NpC>2) + and (Pu0 2 ) ++ we were forced to 
reduce this value substantially in order to minimize 
the rms deviations. It is difficult to estimate the 
accuracy of the determination of Va because the as- 
sumption of sixfold symmetry about the axis may be 
a drastic simplification. 

Similarly we started the calculation with Z ~ 35,000 
cm -1 and were forced to raise it to mimimize the rms 
deviations. Of all the levels below 25,000 cm -1 
i( t>2 (2,1) depends most sensitively on Z. However 
A£ ( 3 4>2) ~ 0.03 AZ so one sees that Z is probably 
not very closely determined by the positions of these 
low levels. 

The values given for X and Y depend on whether 
we have correctly located and identified the broad, 
strong peaks in the absorption spectra. If we assume 
that both the locations and the identifications are 
correct then the accuracy of determination of X and Y 
is probably ±300 cm -1 . 

For the (PuC>2) + + ion the rms deviation between the 
calculated and the observed positions of 12 peaks is 
150 cm -1 . The deviation was calculated from the 
formula: 



A 2 = 2(A£,f/(«-8) 



n is the number of observed levels and 8 is the number 
of parameters. Actually this formula may exaggerate 
the deviation because physical requirements limit 
one's freedom in adjusting the parameters. Also, 
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a more refined variational calculation than we used 
would undoubtedly enable one to choose parameters 
that would decrease the rms deviation somewhat. 
The 150 cm -1 deviation is comparable with the devia- 
tion of 125 cm -1 that Sugar [22] gives for the levels 
of the free Pr 3+ ion. 

For the (Np0 2 ) + ion the rms deviation between the 
calculated and the observed positions of the 14 peaks 
used to fix the parameters is 272 cm -1 . If we include 
also the 3 A3, ^K (2, - 1) and ;i iTo (2, - 1) peaks the rms 
deviation is 241 cm 1 . 

6. Conclusion 

We have shown how the solution absorption spec- 
tra of the (Np02) + and (PuC>2) ++ ions can be inter- 
preted on the basis of exact calculations and plausible 
arguments concerning the strength and width of the 
absorption peaks. To confirm the interpretation it 
would be of some value to locate experimentally the 
two peaks for (NpC>2) + that lie below 6,000 cm -1 and 
the three peaks for (Pu0 2 ) ++ that lie below 9,000 
cm -1 . Further to check the theory developed here 
it is highly desirable to obtain the low temperature 
absorption spectra of a crystalline plutonyl compound 
and of a Np (V) compound. A study of the Zeeman 
effect would be helpful. 

We think the parameters and their variations among 
the ions (Np0 2 ) ++ , (Np0 2 ) + and (Pu0 2 )++ are now 
sufficiently well known so that a theoretical investiga- 
tion of the absorption spectrum of the (Am0 2 ) ++ ion 
could be undertaken. Also, by combining the infor- 
mation gained from studying NpF 6 and the (Pu0 2 ) ++ 
ion, it should be possible to interpret the spectrum 
of PuF 6 . 



We are indebted to Mr. Walter Lipton who coded 
the matrices for the machine calculations, and to Dr. 
W. C. Waggener who sent us a large scale plot of his 
neptunium data. 

7. Appendix 

In analyzing the spectra it is useful to have some 
quantitative ideas, however crude they may be, con- 
cerning the relative breadths and strengths of the 
absorption peaks that correspond to the various 
calculated energy levels. One then has some basis 
for judging whether or not it is reasonable to associate 
a calculated level with an observed absorption peak. 
We have developed a method which gives us some 
insight into this problem. 

In the ground state the electrons are predominantly 
in states with mi =± 3 and ±2. If in the excited state 
the electrons are still predominantly in the ±3 and 
±2 states we assume that the absorption peak will 
be narrow. On the other hand, if in the excited state 
there is an appreciable probability that the mi = ± 1, 
orbitals are occupied, then we assume the absorption 
peak will be broad. We can calculate from the 



eigenvectors the distribution of charge over the 
various orbital states. These calculations have been 
made and the results are given in table 2 for the 
(Pu0 2 ) + + ion and in table 3 for the (Np0 2 ) + ion. 

Now suppose we have two states M* and 4> between 
which a transition can take place. Due to the various 
interactions each one of these states will be a compli- 
cated mixture of the zero order states. Let us there- 
fore write 

*P — a\\f a + bipb + cijjc + . . . . 

<£ = ai// a + /3i///3 + yi//y+. . . . 

Here the i// a , . . ., i// a , . . . denote the zero order 
states and a, . . ., a, . . . are appropriate numerical 
coefficients. Suppose the transition from \jj a to i//« is 
"allowed"; that is, suppose t// rt and i// a differ in, at 
most, the value of mi for one electron. For example, 
i//„ might be 3 H 4 |3-, 2-> and i// a might be 3 n 2 |-3, 2->. 
Similarly, suppose the transition from i//& to i//# is 
"allowed," etc. We then take as a quantitative esti- 
mate of the strength of the peak corresponding to 
the *f? to <t> transition 

a 2 a 2 (i// a -* i// a ) + b 2 (3 2 (tfi b -> fo) + . . . . 

The bracketed quantities denote the probabilities for 
the indicated transitions. 

This device is useful in estimating the relative 
strengths of transitions that would be forbidden if the 
states were pure, but which are somewhat "allowed" 
because of the complicated mixing of the states which 
actually occurs. The numbers a 2 a 2 , 6 2 /3 2 , etc., have 
been calculated and are given in tables 2 and 3 when 
they exceed 0.01. 

The estimates of relative strengths obtained in this 
way can conceivably be improved by using empirical 
information which may be available. For example, 
the strengths of the observed neptunyl lines at 6752 
and 8168 cm -1 , calculated according to the formula 
in the appendix to III, are about 7.7 and 39, respec- 
tively. The former corresponds principally to a 
3^>— 2 transition: the latter to a 3— » — 3 transition. 
The ratios of the strengths of these two transitions 
should therefore be about 1:5. We would expect the 
relative strengths of the 3— > — 2 and 3— » — 3 transi- 
tions for the (Np0 2 ) + and (Pu0 2 ) ^ ions also to be in 
the approximate ratio 1:5. The ratio seems to be 
closer to 1:3 but the difference is perhaps not sur- 
prising in view of the crudeness of the assumptions 
made. Also there may be real differences between the 
one- and the two-electron cases. 
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